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TWO-DIMENSIONAL  FCT  MODEL  OF 
LOW -ALTITUDE  NUCLEAR  EFFECTS 

I.  INTRODUCTION 

Low-yield  experiments  with  chemical  explosives  detonated  at  low  altitude 
have  revealed  unusual  Mach  stem  structures.  Existing  codes  employed  to 
simulate  these  phenomena  have  been  unable  to  resolve  them  satisfactorily. 
During  the  year  beginning  1  June  1979,  personnel  of  the  Laboratory  for 
Computational  Physics  have  been  working  under  a  contract  with  the  Defense 
Nuclear  Agency  to  develop  software  for  performing  these  calculations.  For 
this  purpose  we  have  employed  the  method  of  Flux-Corrected  Transport  (FCT) , 
developed  at  NRL  by  J.  P.  Boris  and  D.  L.  Book. 

The  thrust  of  the  program  during  the  first  year  has  been  toward  testing 
the  various  alternative  forms  of  FCT  which  have  been  developed,  with  a  view 
to  designing  an  optimum  code  (or  codes)  for  solving  alrblast  problems.  The 
principal  features  we  have  tried  to  investigate  and  evaluate  are  (1)  adaptive 
regriddlng;  (2)  multidimensional  (i.e.,  not  time-split)  flux  limiting;'* 

(3)  the  use  of  ZIP  differencing  in  calculating  fluxes,  (4)  the  use  of 
spatially  higher-order  (6th-order,  8th-order,  etc.)  difference  schemes  as  a 
basis  for  designing  the  FCT  routine.  We  concluded  that  (1)  and  (2)  are 
definitely  efficacious;  (1)  indeed  is  indispensable  for  some  problems.  We 
found  strong  indications  that  (3)  is  ineffective  on  the  problems  studied, 
while  we  have  thus  far  not  tested  (4)  and  have  only  indirect  evidence  in 
favor  of  using  high-order  methods.  In  this  report  we  will  discuss  only  cal¬ 
culations  related  to  testing  (1)  and  (2) . 

The  plan  for  FY-80  listed  four  milestones:  (1)  Implement  adaptive  rezone 

in  Cartesian  coordinates;  (2)  Apply  code  to  model  test  problem  (plane  shock 

Manuscript  submitted  August  26,  1980 


1 


above  ramp);  (3)  Convert  code  to  cylindrical  coordinates;  (4)  Apply  code  to 
model  test  problem  (0.5  ton  HE  at  5  m  height). 

All  of  these  have  been  carried  out  and  are  described  below.  In  addition, 
we  have  implemented  a  vectorized  real-air  equation-of-state  routine,  permitting 
us  to  estimate  the  importance  of  real  gas  effects  on  alrblast  calculations. 

Section  II  of  this  report  describes  the  code  development  work  we  have 
performed  [items  (1)  and  (3)  in  the  above  list].  Section  III  describes  our 
shock-on-wedge  calculations  [item  (2)1.  Section  IV  describes  our  work  on 
height-of -burst  studies;  however,  the  parameters  employed  in  the  test  cases 
differ  from  those  specified  above.  The  text  of  a  paper/  describing  some  of 
the  results  obtained  in  this  program,  presented  at  the  Seventh  International 
Conference  on  Numerical  Methods  in  Hydrodynamics  held  during  25-27  June  in 
Stanford,  is  Included  as  an  appendix  to  the  report. 

Work  is  continuing  during  the  remainder  of  the  fiscal  year  to  extend 
and  refine  the  results  reported  here. 


II.  CODE  MODIFICATION  AND  DEVELOPMENT 


A.  Leapfrog  —  Trapezoidal  Code 

A  2D  Cartesian  leapfrog-trapezoidal  code  written  under  DNA  support  by 
Zalesak  for  ionospheric  striation  studies  was  used  as  the  starting  point  in 
developing  a  reflecting  shock  code.  This  was  done  by  disabling  certain 
features  (magnetic  fields,  gravitational  potential)  specific  to  the  original 
problem,  and  writing  initialization  and  boundary  condition  routines  appro¬ 
priate  to  the  shock  problem. 

In  its  present  form  the  code  solves  the  ideal  fluid  equations  on  a 
fixed  uniform  mesh,  using  an  unsplit  2D  FCT  transport  routine.  This  routine 
can  embody  either  a  leapfrog  or  a  trapezoidal  finite-difference  scheme.  The  former 
is  (linearly)  reversible  and  time-centered  and  has  unit  linear  amplification,  but  is 
subject  to  weak  (grid-separation)  instability.  This  is  corrected  by  calling 
the  trapezoidal  scheme  at  intervals  (typically  every  five  timesteps)  to 
"merge"  the  separating  solutions. 

A  perfectly  reflecting  boundary  condition  applied  at  the  bottom  of 

the  mesh  represents  the  wedge  surface;  the  incident  shock  front  makes  an  angle 

0  with  the  vertical,  equal  to  the  wedge  angle.  On  the  other  boundary  cells, 
w 

values  of  the  density  P,  pressure  p,  and  velocity  v  are  assigned  their  values 
in  the  ambient  gas  ahead  or  behind  the  incident  shock,  depending  on 
where  they  are  located  with  respect  to  the  known  shock  position  at  that 
time. 

Results  of  the  calculations  (arrays  of  the  fluid  quantities)  are  dumped 
on  magnetic  tape  at  intervals  of  ten  timesteps.  Graphics  (contour  plots,  time 
histories  of  sensors  at  fixed  locations,  etc.)  are  obtained  by  postprocessing 
the  data  tape. 
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B.  FAST2D 


FAST2D  is  a  two-dimensional  timesplit  FCT  fluid  code  developed  by 

8 

Boris  for  reactive  flow  and  hydrodynamic  stability  calculations.  For 
the  present  task  it  was  utilized  in  a  configuration  very  similar  to  that  of 
the  Leapfrog- Trapezoidal  code.  Initial  and  boundary  conditions  were 
applied  in  identical  fashion. 

The  FAST2D  chemistry  routines  were  omitted,  but  the  dynamic  rezone 
capability  was  retained  in  most  of  the  calculations.  The  latter  exploits 
the  improved  resolution  obtainable  by  f ine-gridding  the  portion  of  the 
mesh  containing  the  most  Important  or  most  intransigent  part  of  the  calcu¬ 
lation  — in  this  case,  the  region  in  which  the  Mach  stem(s)  and  triple 
point(s)  are  located  (Fig.  1).  Coarse  gridding,  which  reduces  the  core 
and  running-time  requirements,  provides  adequate  resolution  on  the  less 
sensitive  remainder  of  the  mesh. 

The  most  important  feature  of  the  rezone  is  its  ability  to  track  the 
critical  region.  This  means  that  the  calculation  is  effectively  carried  out 
in  the  reference  frame  of  the  leading  Mach  stem.  In  calculations  with  glancing 
incident  shocks,  where  the  opening  angle  a  is  as  small  as  one  or  two  degrees, 
the  triple  point  must  move  a  long  way  (several  hundred  cells  on  a  uniform  mesh) 
before  it  is  far  enough  above  the  wedge  (5-10  cells)  to  be  resolved  satis¬ 
factorily.  With  the  dynamic  rezone,  only  ~  150  cells  need  to  be  used  in 
the  calculation. 

C.  Real-Air  Equation  of  State 

A  routine  provided  by  A.  Kuhl  (RDA)  which  calculates  the  effective  ratio  of 

specific  heats  (adiabatic  index)  if  In  a  variety  of  situations  of  interest  for  alrblast 
9 

problems  has  been  converted  to  run  on  the  NRL  TI/ASC.  In  its  original  form 
this  routine  accepts  values  of  specific  internal  energy  e  and  density  P 
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and  returns  y,  temperature  T  and  pressure  p  according  to  four  prescriptions: 
(1)  a  table  lookup  for  dry  real  air,  valid  from  room  temperature  to  ~  10^K 
and  from  millitorr  to  kilobar  pressures;  (2)  a  polynomial  fit  to  these  data; 
(3)  an  analogous  formula  for  detonation  products  from  the  high  explosive  Atlas 
Aquanel;  and  (4)  the  values  appropriate  to  an  ideal  gas  (constant  gamma). 

The  NRL  version  is  vectorized;  it  accepts  arrays  of  values  of  e  and  P  and 
returns  corresponding  arrays  of  y  »  T,  and  p.  The  resulting  code  runs 
approximately  two  orders  of  magnitude  faster  in  mode  (1)  than  the  original 
version,  and  has  been  used  with  both  the  Leapfrog -Trapezoidal  and  FAST2I) 
codes  to  calculate  the  "real  air"  cases  described  below. 
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III.  SHOCK-ON-WEDGE  CALCULATIONS 

This  section  briefly  describes  the  results  of  the  principal  calculations 
we  have  performed  on  shocks  reflecting  from  wedges.  The  runs  reported  here 
(summarized  in  Table  1)  are  discussed  in  three  categories:  small  wedge  angles, 
large  (i.e.,  >  4^)  wedge  angles,  and  runs  carried  out  for  the  purpose  of 
comparing  two  codes  or  models.  The  types  of  Mach  reflection  to  be  expected 
can  be  found  from  Fig.  2,  which  divides  the  parameter  space  into  seven  distinct 
regions . 

A.  Small  Wedge  Angles 

For  small  wedge  angles  (i.e. ,  glancing  incident  shocks),  the  opening  angle 

a  is  ~  10°  and  resolution  of  the  Mach  shock  region  presents  no  great 

difficulty.  The  Leapfrog-Trapezoidal  code  was  used  to  perform  two  calculations 

of  this  sort  (filled  circles  in  Fig.  2),  both  with  an  ideal-gas  EoS.  The 

wedge  angles  were  9  *  26.56°  and  0  *  20°  ;  the  Mach  numbers  were  M  ■  8.06 

w  w 

and  M  »  6.9,  respectively;  and  the  adiabatic  index  y  was  set  equal  to  1.35, 
roughly  the  average  of  the  real-air  values  over  the  range  of  physical  tempera¬ 
tures  and  pressures  encountered  in  these  calculations. 

Results  are  shown  in  Figs.  3-6.  Figure  3a  displays  calculated  density  con¬ 
tours  after  1000  timesteps  for  an  incident  shock  I  arriving  from  the  right  with 
0w  *  20°,  and  M  ■  6.9,  a  case  for  which  complex  Mach  reflection  results 
(Cf.  Fig.  2).  The  incident  shock  I,  the  Mach  stem  ,  the  slip  or  con¬ 

tact  surface  CS,  and  R  and  R',  the  shocks  reflecting  from  the  wedge  and  the 
first  Mach  stem,  respectively,  are  marked.  The  fanlike  lines  labeled  M2  have 
been  added  where  the  density  is  steepening  up,  suggesting  the  system  is  close  to 
developing  a  second  Mach  shock.  The  contour  levels  are  identical  with  those  chosen 
by  Ben-Dor  ^  in  processing  his  interferogram  data  (Fig.  4a). 
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Figure  3b  displays  calculated  density  contours  after  800  timesteps  for  the 

case  with  Q  =  26.56°  and  M  *  8.06.  The  structure  obtained  is  character- 
w 

istic  of  double  Mach  reflection.  The  incident  shock  I,  the  Mach  stems  and 
M2»  the  contact  surface  CS,  and  R  and  R'  are  marked. 

By  comparing  Figs.  3a  and  4a,  we  see  that  not  only  the  shape  of  the 
reflected  wave  fronts,  but  also  the  calculated  contour  levels  interior  to 
them  are  in  excellent  quantitative  agreement  with  their  measured  counterparts. 
Figure  4b  displays  the  calculated  density  and  pressure  profiles  along  the  wedge 
for  this  case,  showing  the  steep  jump  at  the  position  of  the  Mach  stem  M^  and 
the  locations  of  the  two  maxima. 

In  Fig.  5  the  pressure  calculated  on  the  wedge  is  plotted  and  compared  with 
Ben-Dor's  10  experimental  data.  The  discrepancies  are  concentrated  in  the 
two  pressure  peaks,  both  of  which  lie  substantially  above  the  measurements. 

There  are  two  explanations  for  this:  (1)  at  high  temperatures,  the  effective 
value  of  y  in  real  air  drops  (see  Section  IIIC) ,  reducing  the  pressure  jump 
across  a  shock;  and  (2)  Ben-Dor's  data  were  apparently  less  well  resolved 
than  the  numerical  results  (or  possibly  the  peaks  were  eroded  in  the  course 
of  his  data  reduction). 

Figure  6  shows  an  expanded  view  of  the  flow  field  in  the  neighborhood 
of  the  region  where  the  contact  surface  approaches  the  wedge.  Each  line 
segment  denotes  the  magnitude  and  direction  of  the  velocity  defined  at  the 
point  it  is  directed  away  from.  The  velocities  are  viewed  in  the  reference 
frame  moving  with  the  first  Mach  stem;  M2  and  CS  have  the  same  meanings  as 
before.  The  slip  line  has  been  extended  by  following  the  flow  markers.  The 
"jelly-roll"  effect  is  very  pronounced;  two  complete  circuits  can  be 


distinguished.  The  ability  to  resolve  such  structures  is  largely  respon¬ 
sible  for  the  relative  efficacy  of  FCT  algorithms  in  solving  this  kind  of 
problem. 

All  of  the  calculations  on  reflections  from  small-angle  wedges  carried 

out  with  the  Leapfrog-Trapezoidal  code  utilized  uniform  mesh  spacing  with 

fbc  ■  fiy  .  The  meshes  used  extended  up  to  150  x  30  in  size.  Typical  running 
3 

times  for  10  cycles  were  ~  270  s,  i.e.,  ~  60  pa  per  cell-cycle. 

B.  Large  Wedge  Angles 

As  the  wedge  angle  increases  (i.e.,  as  the  incident  shock  strikes  the 
reflecting  surface  more  nearly  normally),  the  opening  angle  a  decreases. 
When  it  goes  to  zero,  we  have  regular  reflection  (region  RR  of  Fig.  2). 
However,  if  a  is  small  but  finite,  then,  as  was  explained  in  Section  IIB, 
resolution  of  the  Mach  reflection  region  is  impractical  unless  a  regridding 
algorithm  is  employed.  The  dynamic  regridding  procedure  of  FAST2D  was  used 
for  this  purpose  in  two  calculations:  0w  =  44°,  M  *  5.074,  and  =  46.5°, 
M  *  5.046,  both  with  the  real  air  EoS.  These  cases  are  indicated  with  open 
circles  in  the  double  Mach  reflection  region  of  Fig.  2.  Note  that,  in 
contrast  with  the  small-wedge-angle  cases,  they  correspond  to  detached 
reflection.  Our  boundary  conditions  are  inconsistent  with  detachment,  but 
this  shduld  not  affect  the  veracity  of  the  calculated  solution  away  from  the 
corner  of  the  wedge. 

Figure  7  shows  pressure  contours  from  four  different  times  in  the  44° 


calculations.  We  have  used  enough  contours  to  show  the  coarse  structure  but 


require  a  larger  number  to  resolve  the  feature  corresponding  to  CS  in  Fig.  3. 
Only  the  finely-gridded  portion  of  the  mesh  is  shown  in  the  figure.  Little 
detail  can  be  discerned  in  frames  (a)  and  (b) .  By  frame  (c) ,  M^  is  distinguish¬ 
able;  in  frame  (d) ,  R,  R',  M^,  and  M2  are  all  well  resolved.  The  contours  in 
the  first  three  frames,  which  were  generated  using  an  ideal-gas  EoS,  are 
sensibly  identical  with  those  at  the  corresponding  times  in  the  real-air  EoS 
calculation  used  to  obtain  the  last  frame. 

Figure  8  shows  plots  of  the  pressure  loading  on  the  wedge  in  orthographic 
projection  at  various  times.  The  contours  tend  to  sharpen  up  at  later  times 
as  resolution  improves;  to  this  extent  the  expected  self-similarity  in  the 
solution  fails  to  materialize.  In  the  final  frame  the  leading  Mach  stem  is 
about  10  cells  high;  to  reach  this  stage  it  had  to  traverse  the  equivalent 
of  450  uniform  grid  points.  Similar  contours  were  produced  in  the  46.5° 
calculations . 

Because  the  0  *  44°  and  0  =46.5°  cases  are  so  similar,  we  discuss 

w  w 

only  the  latter.  The  computed  value  of  a  for  0^  =  46.5°  is  approximately 
2.5°  for  a  real-air  EoS  (with  V  held  constant  at  the  value  1.35  we  obtain 
a»  3.2°).  Either  value  agrees  with  that  measured  to  within  the  experimental 
error, ~  2°.  In  Fig.  9,  we  compare  the  calculated  (using  the  real-air  EoS) 
and  measured  values  of  the  pressure  at  the  surface  of  the  wedge.  The  upper 
traces  are  pressure  plots  vs.  positions  at  two  times  (the  unit  of  time  is 
arbitrarily  designated  to  be  ys ,  whereupon  the  distances  in  cm  become  deter¬ 
mined)  .  The  lower  traces  are  time  histories  at  two  stations  used  in  the  BRL 
experiments  ^  .  (If  the  wave  evolution  adhered  exactly  to  the  approximate 

similarity  description,  spatial  profiles  and  time  histories  would  be  completely 
equivalent) . 
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The  agreement  in  the  shapes  of  the  profile  is  striking.  The  values  of 
the  lower  pressure  peak,  corresponding  to  the  Mach  stem  M^,  are  nearly  identical. 
The  calculated  value  of  the  second  pressure  peak  is  11%  lower  than  the  measured 
value,  and  thus  is  still  within  experimental  uncertainty.  The  rise  in  the 
values  of  both  peaks  and  the  sharpening  in  detail  from  the  first  to  the  second 
plot  probably  has  the  same  explanation  in  both  cases,  namely,  improved  resolu¬ 
tion  as  the  mach  stem  grows.  That  is,  we  believe  that  the  experimental  histories 
are  at  the  limit  of  the  resolution  Imposed  by  the  finite  size  and  time  response 
of  the  pressure  sensors,  which  play  the  same  role  as  space  and  time  truncation 
errors  in  the  finite-difference  approximation.  Figure  10  shows  the  history 
of  the  peak  pressure  on  the  wedge,  this  time  calculated  for  an  ideal  gas  with 
X  =  1.35.  The  effect  of  improvement  in  resolution  with  increasing  time  is 
again  apparent. 

Figure  11  shows  how  well  the  double  Mach  reflection  structure  is  resolved 
in  the  simulation.  The  pressure  in  bar  is  plotted  along  the  trajectory  pressure 
discontinuity  (I,  R'’,  M2,  and  the  flow  "stagnation"  point  on  the  boundary)  in 
succession  transversely.  As  can  be  seen,  the  pressure  "discontinuities"  are 
resolved  equally  well,  within  about  2  cells. 

Movies  have  been  made  of  the  evolving  density  and  pressure  contours  for  the 
0w  =  26.56°  and  =  46.5°  cases, respectively,  the  various  levels  being 
identified  by  different  colors.  These  movies  (several  frames  from  which  are 
reproduced  in  Figs.  7  and  8)  show  graphically  how  the  various  features  of  the 
reflection  structure  sharpen  up  as  sufficient  resolution  becomes  available. 

They  also  show  how  efficiently  the  regridding  algorithm  moves  to  keep  the  re¬ 
flection  region  in  the  finely  gridded  part  of  the  mesh,  tracking  it  smoothly 
without  sudden  or  wide  shifts. 

The  46.5°  movie  demonstrates  the  slow  development  of  the  Mach  structure 
and  permits  us  to  trace  the  locations  of  the  pressure  peaks  along  the  wedge 
surface.  The  26.56°  movie  uses  Ben-Dor’s  density  contours  and  clearly  resolves 
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Che  Double  Mach  structure  by  permitting  the  structure  to  move  twice  as  far 
along  the  wedge  as  in  the  above  test  calculations.  This  capability  results 
directly  from  our  use  of  t.he  adaptive  regrldding  algorithm. 

C.  Code  Comparisons 

A  number  of  trial  calculations  were  performed  during  the  process  of 

code  development  in  order  to  test  the  boundary  and  initial  conditions,  EoS 

routine,  etc.  More  critical  test  runs  were  carried  out  to  compare  the  Leapfrog- 

Trapezoidal  and  FAST2D  codes  under  identical  conditions,  and  to  determine  the 

effect  of  using  the  real-air  EoS  instead  of  the  ideal-gas  EoS. 

For  the  former  test  we  started  with  the  results  obtained  by  the  Leapfrog- 

Trapezoidal  code  for  a  Mach  number  M  =  8.06  and  wedge  angle  0  =  25.56°,  using 

w 

an  ideal-gas  EoS  with  V *  1.35.  Then  we  ran  the  same  problem  using  FAST2D 

with  a  uniform  mesh,  disabling  the  regridding  procedure  and  setting  the  timestep 

equal  to  those  employed  in  the  earlier  run.  Figure  12  shows  the  resulting 

10 

density  contours,  along  with  those  obtained  experimentally  by  Ben-Dor. 
Comparisons  of  the  results  at  identical  time  levels  up  to  cycle  800  yielded 
the  following  conclusions: 

(1)  The  two  codes  yielded  qualitatively  (and  in  most  respects  quantitatively) 
similar  results. 

(2)  Slip  lines  were  better  resolved  by  the  Leapfrog-Trapezoidal  code. 

(3)  The  pressure  calculated  by  the  latter  at  the  location  of  the  wedge  showed 

a  greater  tendency  to  numerical  fluctuations  (Fig.  2).  The  first  pressure 
peak,  just  behind  the  leading  Mach  stem,  was  slightly  higher  in  this  cal¬ 
culation. 

(A)  An  anomalously  low  (less  than  ambient)  pressure  was  produced  at  one  point 
by  FAST2D  because  of  timesplitting.  This,  however,  had  no  apparent  effect 
on  the  rest  of  the  grid,  nor  did  any  discernible  adverse  consequences  of 
timesplitting  appear  elsewhere. 
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To  determine  the  role  played  by  the  real-air  EoS,  we  repeated  the  FAST2D 

runs  described  earlier  (9  =  44°,  M  *  5.074,  and  0  ■  46.5°,  M  =  5.046)  with 

w  w 

constant  V  “  1.35.  (In  order  to  obtain  the  same  pressure  jumps  across  the 
initial  shock,  the  Mach  numbers  had  to  be  adjusted  to  M  *  5.173  and  M  -  5.148, 
respectively).  Agreement  with  the  BRL  data  was  still  good.  At  1000  steps 
both  pressure  peaks  in  the  46.5°  case  were  too  high  in  the  ideal-gas  calculation 
the  one  at  the  Mach  stem  by  ~10%,  and  that  at  the  intersection  of  the  slip 
surface  with  the  wedge  by  ~11%.  Similar  results  were  observed  in  the  44° 
calculations.  The  penalty  for  using  the  real-air  EoS  in  its  optimized  vector 
form  is  about  a  20%  increase  in  running  time. 
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IV.  RELATED  CALCULATIONS 

Several  additional  calculations  were  performed  which  bear  closely  on 
the  airblast  problem:  a  shock  reflecting  from  a  wedge  in  a  reactive  medium, 
and  preliminary  height-of-burst  (HoB)  calculations. 

A.  Reflecting  Shock  in  a  Reactive  Medium 

Under  a  Navy-supported  program  to  study  combustion  hydrodynamics,  LCP 

has  developed  a  simple  numerical  treatment  of  combustion  processes  based  on 

12 

the  induction  time  hypothesis.  This  model  represents  the  chemistry 

through  a  composite  process,  in  which  reactants  begin  to  combine  into  combustion 
products  only  after  a  finite  "induction"  time  has  elapsed.  The  rate  at  which 
the  energy-releasing  reactions  proceed  depends  upon  a  single  parameter,  the 
induction  time.  This  in  turn  is  a  function  of  the  local  thermodynamic  variables. 

As  part  of  the  LCP  Combustion  Hydrodynamics  project,  we  have  applied  the 
induction  time  hypothesis  model  to  study  shock  reflections  in  reactive  gases 
(stoichiometric  mixtures  of  in  air  or  oxygen)  at  low  pressure  (0.1  atm). 

Figure  13  shows  the  time  development  of  a  detonation  initiated  by  a  weak 
reflecting  shock  in  an  H^-air  mixture.  The  sequence  of  six  pressure  contour 
plots  traces  the  evolution  of  a  detonation  wave  initiated  by  complex  Mach 
reflection  at  the  surface  of  the  wedge.  The  incident  shock  was  chosen  so 
that  the  pressure  behind  it  is  too  low  to  cause  detonation  to  take  place 
within  the  time  of  the  calculation.  As  with  the  calculations  of  Section  III, 
we  have  used  open  boundary  conditions  at  the  sides  and  top  of  the  system. 
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Figure  13a  shows  the  pressure  contours  corresponding  to  an  incident 
shock  with  0W  “  25°  and  M  -  4.0  (giving  rise  to  complex  Mach  reflection) 
which  has  just  begun  to  reflect.  The  Mach  stem  is  initially  too  small  to  be 
resolved.  in  frame  b,  the  Mach  stem  becomes  discernible,  but  as  yet  no 
apparent  reaction  has  occurred.  By  frame  (c)  the  material  has  begun  to 
ignite  at  a  position  well  behind  the  current  location  of  the  Mach  stem. 

When  the  Mach  stem  passed  that  position  the  pressure  increase  heated  the 
mixture  sufficiently  to  cause  ignition  after  a  short  induction  time  char¬ 
acteristic  of  the  H2  —  air  mixture.  In  frames  (d)  and  (e) ,  at  later  times, 
the  pressure  at  the  Mach  stem  continues  to  grow,  leading  to  shorter  char¬ 
acteristic  induction  times  for  material  between  the  Mach  stem  and  the  original 
ignition  point.  Thus,  we  see  the  ignited  region  accelerate  along  the  wedge 
surface  toward  the  Mach  stem.  Because  more  energy  is  being  released  as  the 
burning  continues,  the  boundary  of  the  ignited  region  also  accelerates  in 
the  direction  of  the  reflected  shock  and  Mach  stem,  as  we  see  from  the 
decrease  in  the  separation  of  the  pressure  contours  near  both  locations. 

A  stable  detonation  pattern  has  not  yet  emerged,  however.  This  is  evident 
from  the  bending  out  of  the  Mach  stem  and  the  lower  density  of  contours 
between  the  Mach  stem  and  the  reflected  shock/detonation  front. 

A  similar  calculation  has  been  carried  out  in  which  all  four  boundaries 
were  made  reflecting.  This  is  seen  to  be  subcritical,  i.e.,  a  stable  detonation 
wave  does  not  form.  When  the  incident  shock  Mach  number  is  increased  to  M  =  5.0, 
however,  a  stable  detonation  wave  does  occur.  As  it  evolves,  an  intricate 
pattern  of  multiple  reflections  from  the  boundaries  develops.  Color  movies 
have  been  made  of  the  pressure  contours  for  both  calculations,  illustrating 
clearly  the  difference  in  the  two  physical  situations. 


These  studies  have  led  to  the  discovery  of  previously  unreported  pockets 


of  unreacted  gas  which  become  folded  in  behind  the  oscillating  cellular 
multidimensional  structure  which  characterizes  most  real  detonations.  These 
pockets,  which  react  eventually  some  time  after  the  main  shock  and  reaction 
zone  have  passed,  ought  to  be  experimentally  observable.  They  also  may  prove 
to  be  the  mechanism  by  which  detonations  sustain  themselves  at  abrupt  area 
changes  and  in  free  space  where  there  are  no  walls  to  provide  new  cellular 
detonation  centers. 

B.  Height-of -Burst  Calculations 

FAST2D  is  written  in  such  a  way  as  to  permit  switching  from  Cartesian 
to  cylindrical  geometry  by  changing  a  single  input  parameter.  A  sample  HoB 
calculation  was  carried  out  by  introducing  a  1  kg  explosive  charge  at  a  height  0f 
7.5  cm  approximately  reproducing  the  conditions  of  the  Mighty  Mach  experiments. 

The  resulting  blast  wave  is  initially  spherical  (i.e.,  semicircular  when  pro¬ 
jected  on  the  r-£  half  plane;  see  Fig,  i4).By  using  vertically  stretched  cells 
at  the  top  of  the  mesh,  it  is  possible  to  prevent  the  blast  wave  from  reflecting 
back  into  the  ground  reflection  region  until  late  in  the  development  of  the 

expanding  shocked  system. 

Unlike  the  constant-strength  incident  shocks  discussed  previously,  the 

pressure  (or  density,  etc.)  profile  through  the  front  varies  temporally  in 

two  ways:  it  weakens  due  to  spherical  expansion,  and  it  steepens  and  narrows 

because  the  subsequent  rarefaction  wave  tends  to  catch  up  with  the  shock. 

A  third  effect  contributes  significantly  to  time  variation  in  real  airblasts 
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produced  by  high  explosives:  the  delayed  combustion  of  the  HE  charge.  This 

could  be  properly  treated  using  a  simple  model  like  that  based  on  the  induction 
induction  time  hypothesis  (Section  IV  A  ) ,  and  a  turbulent  mixing  model,  coupled 
with  an  accurate  EoS  for  the  detonation  products.  The  present  calculation  neglected 
the  effect  of  the  delayed  energy  release, which  will  be  included  in  the  future. 


Because  the  incident  wave  structure  differs  in  both  space  and  time  from 
that  assumed  in  the  shock-on-wedge  problems,  the  reflected  wave  now  looks 
very  different  (Fig.  15).  The  Mach  stem  is  still  visible,  but  the  rare¬ 
faction  wave  prevents  the  formation  of  a  second  Mach  shock  and  distorts 
the  reflected  wave  and  contact  surface. 

Finally,  as  an  example  of  what  can  be  done  with  the  computational  tools 
in  hand,  we  performed  a  HoB  calculation  with  a  ground  structure.  The  latter 
was  a  block  having  a  front  face  sloping  at  45°  (Fig.  14).  Reflecting 
boundary  conditions  were  employed  on  the  structure,  with  the  ramp  modeled  by 
a  stair-step  boundary.  In  addition  to  contour  plots  like  those  previously 
described,  pressure  histories  were  obtained  at  four  locations  along  the 
ramp  (Fig.  16).  The  traces  on  the  left  in  Fig.  16  represent  time  histories 
of  the  maximum  (as  a  function  of  horizontal  position)  pressures  on  the  mesh 
at  the  elevations  corresponding  to  sensors  1,  2,  3,  4,  respectively,  showing 
that  these  increase  with  elevation.  It  is  interesting  to  note  that  the  highest 
peak  pressure  on  the  structure  was  recorded  at  location  2,  not  at  the  bottom 


of  the  ramp. 
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COARSE  ZONES  30  FINE  ZONES 


schematically  are  incident,  reflected  and  Mach  shocks  (solid  lines)  and  slip  surface 
(broken  lines)  for  incident  shock  arriving  from  the  left.  The  finely  gridded  region 
moves  with  Mach  stems  and  triple  points  to  provide  resolution  where  it  is  needed  mos 


WEDGE  ANGLE  (DEGREES) 


MACH  NO.  - »- 

Fig.  2  -  Types  of  shock  reflections:  RR,  SMR,  CMR,  and 
DMR  denote  regular  reflection  and  single,  complex  and  dou¬ 
ble  Mach  reflection,  respectively.  D  and  A  label  detached 
and  attached  shocks,  respectively.  Closed  circles  repre¬ 
sent  calculations  done  with  the  Leapfrog -Trapazoidal  code; 
open  circles  represent  those  done  with  FAST2D. 
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mg.  3  -  Wave  structure  and  density  contours  with  lines  of  discontinuity  identified  for 
(a)  complex  Mach  reflection  case  (9  =  20°,  M  =  6.9);  (b)  double  Mach  reflection  case 

(0  =  26.56°,  M  =  8.06).  W 


REDUCED  DENSITY  CONTOURS  (M,  =  6.9, 0W  =  20°, y  =  1.35)  393.85 
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Fig.  4  -  (a)  Contours  of  Figure  3(a) ,  labeled  with  density  values  used  by  Ben- 
Dor;  (b)  Density  (dots)  and  pressure  (plus  sign)  contours  calculated  on  wedge 
(circles,  fitted  by  solid  line)  and  calculated  value  (dots). 
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Fig.  5  -  Comparison  in  complex  Mach  reflection  case  of  experi¬ 
mentally  measured  density  on  wedge  (circles,  fitted  by  solid 
line)  and  calculated  value  (dots). 
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Fig.  6  -  Rollup  of  slip  line  at  wedge  surface  (solid  line).  Short  line  segments  rep 
resent  magnitude  and  direction  of  velocity  at  mexh  points  (dots),  in  reference  frame 
of  Mach  stem. 
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Pig.  7  -  Development  of  Mach  shock  in  ideal-gas  44-degree_galcula- 
tion.  Pressure  contours  are  shown  at  times_Ca)  1.57  x  10  ;  (b) 

4.09  X  10  ;  (c)  5.95  x  10”°;  (d)  1.54  x  10  (this  last  frame  ta¬ 

ken  from  calculation  using  real  air  EOS). 


Fig.  8  -  Pressure  plots  from  ideal -gas  44 -degree  calcu¬ 
lation  in  orthographic  projection  at  times  (a)  0.00;  (b) 
1.13  x  10”4;  (c)  1.60  x  10-4;  (taken  from  calculation 
using  real-air  EOS) . 
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Fig.  9  -  Upper  and  lower  diagrams  show  pressure  in  PSIA  on  the  wedge  as 
a  function  of  position  fore  =  46.5°  (real-air  EOS)  at  two  times  in  the  simu¬ 
lation  (upper  traces  and  atwtwo  stations  in  tne  BLR  experiments  11(low- 
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Fig.  10  -  Pressure  profile  on  wedge  in  46 . 5° ideal-gas  calculation  at  four  differ 
ent  times,  showing  the  effect  of  improving  resolution. 
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Fig.  11  -  Plot  of  pressure  vs  position  (46.5  wedge  angle)  along 
the  trajectory  1 -2-3-4  shown  in  the  inset  which  intersects  four 
surfaces  of  discontinuity  (solid  lines)  transversely. 


FAST2D 

Fig.  12  -  Comparison  between  experimental,  Leapfrog 
trapezoidal,  and  FAST2D  double-Mach -ref lection  den¬ 
sity  contours. 


13  -  Pressure  contours  produced  in  detonation  calculation  (reflected  shock  in  reac 
medium)  at  times  (a)  1.10  x  lO"6^;  (b)  2.70  x  10“6ns;  (c)  3.62  x  10"6  ns;  (d)  3.91 
is;  (e)  4.50  x  10-6^;  (f)  6.73  x  10“6 


Fig.  14  -  Schematic  of  HoB  calculation.  Energy  release  is  at  dot  on  vertical  line 
(z  axis).  Shocks  (dashed  lines)  reflect  off  boundaries,  limiting  the  extent  of 
validity  in  space-time  of  calculation.  Hatched  region  indicates  location  of  solid 
(reflecting)  ground  structure. 


Fig.  15  -  Development  of  Mach  shock 
In  HoB  calculation.  Pressure  con¬ 
tours  are  shown  at  times  (a)  5.70  x 
10-5;  (b)  7.03  x  IQ"5;  (c)  1.62  x  10”4 


t  (arbitrary  units)  - ►- 

Pressure  histories  at  four  sensor  locations  on  ramp  of  structure  sketched 
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04 cy  Attn  TITL 
Olcy  Attn  DDST 

Defense  Technical  Information  Center 
Cameron  Station 
Alexandria,  VA  22314 

(12  If  open  pub,  otherwise  2  - 
no  WN INTEL) 

12cy  Attn  DD 

Chairmen 

Department  of  Defense  Explo  Safety  Board 
Rm  856-C 

Hoffman  Building  1 
2461  Elsenhower  Avenue 
Alexandria,  VA  22331 
Olcy  Attn  Chairman 

Commander 
Field  Command 
Defense  Nuclear  Agency 
Klrtland  AFB,  NM  87115 
Olcy  Attn  FCTMOF 
Olcy  Attn  FCT 
Olcy  Attn  FCPR 


Chief 

Field  Command 
Defenae  Nuclear  Agency 
Livermore  Division 
F  0  Box  808  L-317 
Livermore,  CA  94550 
Olcy  Attn  FCPRL 

Director 

Joint  Strat  TGT  Planning  Staff 
Offutt  AFB 
Omaha,  NB  68113 
Olcy  Attn  DOXT 
Olcy  Attn  JLA 
Olcy  Attn  JLTW-2 
Olcy  Attn  NRI-STINFO  Library 
Olcy  Attn  EPFS 

Commandant 

Nato  School  (Shape) 

AFO  New  York,  NY  09172 

Olcy  Attn  U.S.  Documents 
Officer 

Under  Secy  of  Def  for  Rsch  & 

Engrg 

Department  of  Defense 
Washington,  DC  20301 

Olcy  Attn  Strategic  & 

Space  Systems  (OS) 

Director 

BMD  Advanced  Technology  Center 
Department  of  the  Army 
F  0  Box  1500 
Huntsville,  AL  35807 
Olcy  Attn  1CRDABH-X 
Olcy  Attn  ATC-T 

Commander 

BMD  Systems  Command 
Department  of  the  Army 
F  0  Box  1500 
Huntsville,  AL  35807 

Olcy  Attn  BMDSC-H  N  HURST 

Chief  of  Engineers 
Department  of  the  Army 
Forrestal  Building 
Washington,  DC  20314 
Olcy  Attn  DAEN-RDM 
Olcy  Attn  DAEN-MCE-D 
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Deputy  of  Chief  of  Staff  for 
OPS  &  Plans 
Department  of  the  Army 
Washington,  DC  20310 
Olcy  Attn  DAMO-NC 

Commander 

Harry  Diamond  Laboratories 
Department  of  the  Army 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 

(CNWDI-Inner  Envelope: 

Attn:  DELHD-RBH) 

Olcy  Attn  Chief  Div  20000 
Olcy  Attn  DELHD-I-TL  (Tech  Lib) 

Commander 

U.S.  Army  Armament  Material  Readiness 
Command 

Rock  Island,  IL  61202 
Olcy  Attn  MA  Library 

Director 

U.S.  Army  Ballistic  Research  Labs 
Aberdeen  Proving  Ground,  MD  21005 
Olcy  Attn  DRDAR-BLE  J  Keefer 
02cy  Attn  DRDAR-TSB-S  (Tech  Lib) 
Olcy  Attn  DRDAR-BLT  W  Taylor 
Olcy  Attn  DRDAR-BLV 

Commander 

U.S.  Army  Communications  Conanand 
Fort  Huachuca,  AZ  85613 

Olcy  Attn  Technical  Reference  Div 

Commander 

U.S.  Army  Concepts  Analysis  Agency 
8120  Woodmont  Avenue 
Bethesda,  MD  20014 

Olcy  Attn  M0CA-ADL  (Tech  Lib) 

Commander 

U.S.  Army  Engineer  Center 
Fort  Belvoir,  VA  22060 
Olcy  Attn  ATZA 

Division  Engineer 
U.S.  Army  Engineer  Div  Huntsville 
P  0  Box  1600,  West  Station 
Huntsville,  AL  35807 
Olcy  Attn  HNDED-SR 

Division  Engineer 

U.S.  Army  Engineer  Div  Ohio  River 

P  0  Box  1159 

Cincinnati,  OH  45201 

(Unclassified  Only) 

Olcy  Attn  0RDAS-L  (Tech  Lib) 


Director 

U.S.  Army  Engr  Waterways  Exper 
Station 
P  0  Box  631 
Vicksburg,  MS  39180 

Olcy  Attn  WESSD  G  Jackson 
Olcy  Attn  WESSA  W  Flathau 
Olcy  Attn  J  Strange 
Olcy  Attn  WESSE  L  Ingram 
Olcy  Attn  Library 

Commander 

U.S.  Army  Foreign  Science  &  Tech  Ctr 
220  7th  Street,  HE 
Charlottesville,  VA  22901 
Olcy  Attn  DRXST-SD 

Commander 

U.S.  Army  Material  &  Mechanics 
Ssch  Ctr 

Watertown,  MA  02172 

(Address  CNWDX:  Attn: 

Document  Control  for) 

Olcy  Attn  Technical  Library 
Olcy  Attn  DRXMR  J  Mescall 
Olcy  Attn  DRXMR-TE  R  SHEA 

Commander 

U.S.  Army  Material  Dev  &  Readiness 
CMD 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Olcy  Attn  DRCDE-D  L  Flynn 
Olcy  Attn  DRXAM-TL  (Tech  Lib) 
Unci  only 

Commander 

U.S.  Army  Missile  Command 
Redstone  Arsenal,  AL  35809 
Olcy  Attn  DRDMI-XS 
Olcy  Attn  RSIC 

Commander 

U.S.  Army  Mobility  Equip  RAD  CMD 
Fort  Belvoir,  MD  22060 

(CNWDI  to  Army  Mat  Dev 
6  Readiness  Conmand) 

Olcy  Attn  DRDME-WC  (Tech  Lib) 

Commander 

U.S.  Army  Nuclear  &  Chemical  Agency 
7500  Backlick  Road 
Building  2073 
Springfield,  VA  22150 

(desires  only  Icy  to  Library) 
Olcy  Attn  Library 


Commandant 

U.S.  Army  War  College 
Carlisle  Barracks,  PA  17013 
Olcy  Attn  Library 

Commander 

David  Taylor  Naval  Ship 
R&D  CTR 

Bethesda,  MD  20084 

(CNWDI  only  Attn 
Mrs.  M.  Blrkhead 
Code  5815.61) 

Olcy  Attn  Code  L42-3 
(Library) 

Of  f lcer-in-Charge 

Naval  Construction  Battalion  Center 
Civil  Engineering  Laboratory 
Port  Hueneme,  CA  93041 

Olcy  Attn  Code  LSI  S  Takahashi 
Olcy  Attn  Code  L51  R  Odello 
Olcy  Attn  Code  L08A  (Library) 

Commander 

Naval  Electronic  Systems  Command 
Washington,  DC  20360 

Olcy  Attn  PME  117-21 

Commander 

Naval  Facilities  Engineering  Command 
Washington,  DC  20390 

Olcy  Attn  Code  09M22C  (Tech  Lib) 
Olcy  Attn  Code  03T 
Olcy  Attn  Code  04B 

Headquarters 

Naval  Material  Comnand 

Washington,  DC  20360 

Olcy  Attn  MAT  08T-22 

Commander 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 

Olcy  Attn  Code  4471  (Tech  Lib) 
Olcy  Attn  Code  013  E  Cooper 

Superintendent 

Naval  Postgraduate  School 

Monterey,  CA  93940 

(desires  no  CNWDI  Documents) 
Olcy  Attn  Code  0142  Library 

Co  man  ding  Officer 
Naval  Research  Laboratory 
Washington,  DC  20375 

(RD  &  RD/N  Attn  Code  1221  for 
&  FRD  Attn  Code  2628  for) 

Olcy  Attn  Code  2627  (Tech  Lib) 
Olcy  Attn  Code  Boris/Book 


Comnander 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

Olcy  Attn  Sea-09G53  (Lib) 

Olcy  Attn  Sea-0351 

Officer  In  Charge 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 
Olcy  Attn  Code  F31 
Olcy  Attn  H.  Gias 

Comnander 

Naval  Surface  Weapons  Center 
Dahlgven,  VA  22448 

Olcy  Attn  Tech  Library  A 
Infor  Services  Brnch 

President 

Naval  War  College 

Newport,  RI  02840 

Olcy  Attn  Code  E-ll 
(Tech  Service) 

Comnander 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Olcy  Attn  Code  266  C  Austin 
Olcy  Attn  Code  3201  P  Cordle 
Olcy  Attn  Code  233  (Tech  Lib) 

Commanding  Officer 
Naval  Weapons  Evaluation  Facility 
Kirtland  Air  Force  Base 
Albuquerque,  NM  87117 
Olcy  Attn  R  Hughes 
Olcy  Attn  Code  10  (Tech  Lib) 

Office  of  Naval  Research 
Arlington,  VA  22217 

Olcy  Attn  Code  474  N  Perron 
Olcy  Attn  Code  715  (Tech  Lib) 
(Unci  only) 

Office  of  the  Chief  of  Naval 
Operations 

Washington,  DC  20350 
Olcy  Attn  OP  981 
Olcy  Attn  OP  03EG 

Director 

Strategic  Systems  Project  Office 
Department  of  the  Navy 
Washington,  DC  20376 

Olcy  Attn  NSP-43  (Tech  Lib) 
Olcy  Attn  NSP-272 


Air  Fores  Geophysics  Laboratory 
Hanscom  AFB,  MA  01731 

Oley  Attn  LWW  K  Thompson 

Air  Force  Institute  of  Technology 
Air  University 

Wrlght-Patterson  AFB,  OH  45433 

(Does  not  desire  classified 
documents) 

Olcy  Attn  Library 

Headquarters 

Air  Force  Systems  Command 
Andrew  AFB 
Washington,  DC  20334 
Olcy  Attn  DLW 

Air  Force  Weapons  Laboratory,  AFSC 
Kir t land  AFB,  NM  87117 

Olcy  Attn  NTES-C  R  Benny 
Olcy  Attn  NT  ED  R  Matalucd 
Olcy  Attn  NTE  M  Flamondon 
Olcy  Attn  R  Guice 
Olcy  Attn  SUL  W  Lee 
Olcy  Attn  DEX 

Assistant  Chief  of  Staff 
Intelligence 

Department  of  the  Air  Force 
Washington,  DC  20330 
Olcy  Attn  IN 

Ballistic  Missle  Office/DE 
Air  Force  Systems  Command 
Norton  AFB,  CA  92409 

(Civil  Engineering) 

Olcy  Attn  DEB 

Ballistic  Missle  Offlce/MN 
Air  Force  Systems  Command 
Norton  AFB,  CA  92409 

(Minuteman)  MNNX 
Olcy  Attn  MNNXH  D  Gage 

Deputy  Chief  of  Staff 
Research,  Development,  &  ACC 
Department  of  the  Air  Force 
Washington,  DC  20330 
Olcy  Attn  AFRDQSM 

Deputy  Chief  of  Staff 
Logistics  &  Engineering 
Department  of  the  Air  Force 
Washington,  DC  20330 
Olcy  Attn  LEEE 


Commander 

Foreign  Technology  Division,  AFSC 
Wrlght-Patterson  AFB,  OH  45433 
Olcy  Attn  NIIS  Library 

Commander 

Rome  Air  Development  Center,  AFSC 
Grlffiss  AFB,  NY  13441 
(Desires  no  CNWDZ) 

Olcy  Attn  TSLD 

Stragetic  Air  Command/XFFS 
Department  of  the  Air  Force 
Offutt  AFB,  NB  68113 

Olcy  Attn  NRI-STINFO  Library 
Olcy  Attn  XPFS 

Department  of  Energy 
Albuquerque  Operations  Office 
F  0  Box  5400 
Albuquerque,  NM  87115 
Olcy  Attn  CTID 

Department  of  Energy 
Washington,  DC  20545 

Olcy  Attn  0MA/RDVE 

Department  of  Energy 
Nevada  Operations  Office 
F  0  Box  14100 
Las  Vegas,  NV  89114 

OJcy  Attn  Mail  &  Records 

for  Technical  Library 

Department  of  the  Interior 
Bureau  of  Mines 
Bldg.  20,  Denver  Federal  Ctr 
Denver,  CO  80225 

Olcy  Attn  Tech  Lib  (Und  only) 

Director 

Federal  Emergency  Management  Agency 
1721  I  Street,  NW 
Washington,  DC  20472 

Olcy  Attn  Hazard  Eval  &  Vul 
Red  Div 

Aerospace  Corp. 

F  0  Box  92957 
Los  Angeles,  CA  90009 
Olcy  Attn  H  Mlrels 
Olcy  Attn  Technical  Infor 
Services 
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Agbablan  Associates 
250  M  Nash  Strait 
El  Sagundo,  CA  90245 

Olcy  Attn  M  Agbablan 

Analytic  Sarvicas,  Inc. 

400  Amy-Navy  Driva 
Arlington,  VA  22202 

Olcy  Attn  G  Hasselbacher 

Applied  Theory,  Inc. 

1010  Westwood  Blvd. 

Los  Angelas,  CA  90024 

(2cys  if  unclass  or 
ley  If  class) 

Olcy  Attn  J  Trullo 

Artec  Associates,  Inc. 

26046  Eden  Landing  Hoad 
Hayward ,  CA  94545 
Olcy  Attn  S  Gill 

Avco  Research  &  Systems  Group 
201  Lowell  Street 
Wilmington,  MA  01887 

Olcy  Attn  Library  A830 

BUM  Corp . 

7915  Jones  Branch  Drive 
McLean,  VA  22102 

Olcy  Attn  A  Lavagnino 
Olcy  Attn  T  Neighbors 
Olcy  Attn  Corporate  Library 

BUM  Corp. 

P  0  Box  9274 
Albuquerque,  KM  87119 
Olcy  Attn  R  Hensley 

Boeing  Co. 

P  0  Box  3707 
Seattle,  WA  98124 

Olcy  Attn  M/S  42/37  R  Carlson 
Olcy  Attn  Aerospace  Library 

California  Research  &  Technology,  Inc. 
6269  Varlel  Avenue 
Woodland  Hills,  CA  91364 
Olcy  Attn  Library 
Olcy  Attn  K  Kreyenhagen 

California  Research  &  Tech,  Inc. 

4049  First  Street 
Livermore,  CA  94550 
Olcy  Attn  D  Orphal 


Calspan  Corp. 

P  0  Box  400 
Buffalo,  NY  14225 
Olcy  Attn  Library 

Denver,  University  of 
Colorado  Seminary 
Denver  Research  Institute 
P  0  Box  10127 
Denver,  CO  80210 

(Only  Icy  of  dess  rpts) 

Olcu  Attn  Sec  Officer  for 
J  Wlsotskl 

EG&G  Washington  Analytical 
Services  Center,  Inc. 

P  0  Box  10218 
Albuquerque,  NM  87114 
Olcy  Attn  Library 

Eric  H.  Wang 

Civil  Engineering  Rsch  Fac 
University  of  New  Mexico 
University  Station 
P  0  Box  25 

Albuquerque,  NM  87131 
Olcy  Attn  N  Baum 

Gard,  Inc. 

7449  N  Natchez  Avenue 
Niles,  IL  60648 

Olcy  Attn  G  Neldhardt 
(Und  only) 

General  Electric,  Co 
Space  Division 
Valley  Forge  Space  Center 
P  O  Box  8555 
Philadelphia,  FA  19101 
Olcy  Attn  M  Bortner 

General  Electric  Co. -Tempo 
816  State  Street  (P  0  Drawer  QQ) 
Santa  Barbara,  CA  93102 
Olcy  Attn  DASIAC 

General  Research  Corp. 

Santa  Barbara  Division 
P  0  Box  6770 
Santa  Barbara,  CA  93111 
Olcy  Attn  B  Alexander 

Higgins,  Auld  Association 
2601  Wyoming  Blvd  NE 
Albuquerque,  NM  87112 
Olcy  Attn  J  Bratton 
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IIT  Research  Institute 
10  W  35th  Street 
Chicago,  IL  60616 

Olcy  Attn  Documents  Library 
Olcy  Attn  R  Welch 
Olcy  Attn  M  Johnson 


Information  Science,  Inc. 
123  W  Padre  Street 
Santa  Barbara,  CA  93105 
Olcy  Attn  W  Dudziak 


Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Olcy  Attn  Classified  Library 


J  H  Wiggins  Co.,  Inc. 

1650  S  Pacific  Coast  Highway 
Redondo  Beach,  CA  90277 
Olcy  Attn  J  Collins 


Raman  Avldyne 
83  Second  Street 
Northwest  Industrial  Park 
Burlington,  MA  01803 
Olcy  Attn  Library 
Olcy  Attn  E  Crlsdone 
Olcy  Attn  N  Hobbs 
Olcy  Attn  R  Ruetenlk 


Raman  Sciences  Corp. 

P  0  Box  7463 

Colorado  Springs,  CO  80933 
Olcy  Attn  F  Shelton 
Olcy  Attn  Library 


Raman  Sciences  Corp. 

Southern  California  Operations 
101  Continental  Blvd  Suite  855 
El  Segundo,  CA  90245 
Olcy  Attn  D  Sachs 


Lockheed  Missiles  &  Space  Co., 
P  0  Box  504 
Sunnyvale,  CA  94086 

Olcy  Attn  TlC-Hbrary 


Inc. 


Los  Alamos  National  Scientific  Lab. 

Mail  Station  5000 

P  0  Box  1663 

Los  Alamos,  NM  87545 

Olcy  Attn  MS  670/J  Hopkins 

Attn  DOC  CON  for  M  Sanpford 
Attn  DOC  CON  for  R  Whittaker 
Attn  DOC  CON  for  MS  364 
(Class  Reports  Lib) 

Attn  DOC  CON  for  6  Spillman 
Attn  DOC  CON  for  A  Davis 
Attn  DOC  CON  for  R  Brldwell 


Olcy 

Olcy 

Olcy 


Olcy 

Olcy 

Olcy 


Lovelace  Biomedical  &  Environmental 
Research  Institute,  Inc. 

P  0  Box  5890 
Albuquerque,  NM  87115 

Olcy  Attn  R  Jones  (Undas  only) 


Martin  Marietta  Corp. 

P  0  Box  5837 
Orlando,  FL  32855 

Olcy  Attn  G  Fotleo 


McDonnell  Douglas  Corp. 

5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
Olcy  Attn  R  Halprln 


Merritt  Cases,  Inc. 

P  0  Box  1206 
Redlands,  CA  92373 

Olcy  Attn  J  Merritt 
Olcy  Attn  Library 


Lawrence  Livermore  Rational  Lab. 

P  O  Box  808 
Livermore,  CA  94550 

Olcy  Attn  DOC  CON  for  L-200  T  Butkovlch 

Olcy  Attn  DOC  CON  for  lech  Inf or  Dept.  Library 

Olcy  Attn  DOC  CON  for  L-205  J  Hearst  (Class  L-203) 

Olcy  Attn  DOC  CON  for  L-90  D  Morris  (Class  L-504) 

Olcy  Attn  DOC  CON  for  L-437  R  Schock 

Olcy  Attn  DOC  CON  for  L-7  J  Hahn 

Olcy  Attn  DOC  CON  for  L-96  L  Woodruff  (Class  L-94) 

Olcy  Attn  DOC  CON  fro  L-90  R  Dong 


Meteorology  Research,  Inc. 

464  W  Woodbury  Road 
Altadena,  CA  91001 
Olcy  Aten  W  Green 

Nathan  M.  Newark  Consult 
Eng  Services 
B106A  Civil  Eng  Bldg. 

University  of  Illinois 
Urbans,  IL  61801 

Olcy  Attn  N  Newark 

Oak  Ridge  National  Lab. 

Nuclear  Division 
Z-10  Lab  Records  Div 
P  0  Box  X 

Oak  Ridge,  TN  37830 

Olcy  Attn  Civil  Def  Res  Proj 
Olcy  Attn  DOC  CON  for  Central 
Research  Library 

Pacifica  Technology 
P  0  Box  148 
Del  Mar,  CA  92014 
Olcy  Attn  G  Kent 
Olcy  Attn  R  Bjork 

Physics  International  Co. 

2700  Merced  Street 
San  Leandro,  CA  94577 
Olcy  Attn  E  Moore 
Olcy  Attn  L  Behrmann 
Olcy  Attn  Technical  Library 
Olcy  Attn  7  Sauer 

R  &  D  Associates 
P  0  Box  9695 

Marina  Del  Rey,  CA  90291 

Olcy  Attn  Technical  Infor  Ctr 

Olcy  Attn  A  Latter 

Olcy  Attn  A  Kuhl 

Olcy  Attn  J  Carpenter 

Olcy  Attn  C  MacDonald 

Olcy  Attn  R  Port 

Olcy  Attn  J  Levis 

Rand  Corp. 

1700  Main  Street 
Santa  Monica,  CA  90406 
Olcy  ATTN  C  Mow 

Sand la  Laboratories 
Livermore  Laboratory 
P  0  Box  969 
Livermore,  CA  94550 

Olcy  Attn  DOC  CON  for  Library 
4  Security  Class  Division 


Sandia  National  Laboratories 
P  0  Box  5800 
Albuquerque,  NM  87185 

(Attn  Mail  Services  Section 
for  Intended  Recipient) 

Olcy  Attn  Mall  Ser  See  W  Roherty 
Olcy  Attn  Mall  Ser  Sec  3141 
Olcy  Attn  Mall  Ser  See  L  Vortman 
Olcy  Attn  Mail  Ser  Sec  A  Chaban 
Olcy  Attn  Mail  Ser  Sec  L  Hill 

Science  Applications,  Inc. 

P  0  Box  2351 
La  Jolla,  CA  92038 

Olcy  Attn  Technical  Library 

Science  Applications,  Inc. 

1257  Tasman  Drive 
Sunnyvale,  CA  94086 
Olcy  Attn  J  Dlshon 

Science  Applications,  Inc. 

2450  Washington  Avenue 
San  Leabdro,  CA  94577 
Olcy  Attn  D  Maxwell 
Olcy  Attn  D  Bernstein 

Science  Applications,  Inc. 

P  0  Box  1303 
McLean,  VA  22102 

Olcy  Attn  M  Knasel 
*  Olcy  Attn  B  Chambers  III 
Olcy  Attn  R  Slevers 
Olcy  Attn  J  Cockayne 

Southwest  Research  Institute 
P  0  Drawer  28510 
San  Antonio,  TX  78284 
Olcy  Attn  W  Baker 
Olcy  Attn  A  Wenzel 

SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

Olcy  Attn  G  Abrahamson 

Systems,  Science  6  Software,  Inc. 

P  0  Box  1620 
U  Jolla,  CA  92038 
Olcy  Attn  Library 
Olcy  Attn  D  Grine 
Olcy  Attn  T  Riney 
Olcy  Attn  K  Pyatt 


41 


I 


Teledyne  Brown  Engineering 
Cunnings  Research  Perk 
Huntsville,  AL  35807 

Oley  Attn  J  Ravenacraft 

Terra  Tek,  Inc. 

420  Uekera  Hey 
Salt  Lake  City,  UT  84108 
Olcy  Attn  Library 
Olcy  Attn  S  Creen 
Olcy  Attn  A  Jones 

Tatra  Tech,  Inc. 

630  N  Rosemead  Blvd. 

Pasadena,  CA  91107 
Olcy  Attn  L  Hwang 
Olcy  Attn  Library 

TRW  Defense  &  Space  Sys  Group 
One  Space  Park 
Redondo  Beach,  CA  90278 
Olcy  Attn  I  Alber 
Olcy  Attn  Tech  Infor  Ctr 
02cy  Attn  N  Lipner 
Olcy  Attn  P  Bhutta 

,  Olcy  Attn  D  Baer 

Olcy  Attn  R  Plebuch 

TRW  Defense  &  Space  Sys  Group 
P  0  Box  1310 

San  Bernardino,  CA  92402 
Olcy  Attn  E  Hong 
Olcy  Attn  P  Dai 

Universal  Analytics,  Inc. 

7740  W  Manchester  Blvd 
Playa  Del  Rey,  CA  90291 
Olcy  Attn  E  Field 

Weidllnger  Assoc.,  Consulting  Eng 
110  E  59th  Street 
Hew  York,  NY  10022 
Olcy  Attn  M  Baron 

I  Weidllnger  Assoc.,  Consulting  Eng 

3000  Sand  Hill  Road 
Menlo  Park,  CA  94025 
Olcy  Attn  J  Isenberg 

Westlnghouse  Electric  Corp. 

Marine  Division 
Handy  Avenue 
Sunnyvale,  CA  94088 
Oley  Attn  W  Vols 


42 


